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Abstract 



1 We describe interferometric observations of the asteroid (41) Daphne in 

2 the thermal infrared obtained with the Mid- Infrared Interferometric Instrument 

3 (MIDI) and the Auxiliary Telescopes (ATs) of the European Southern Obser- 

4 vatory (ESO) Very Large Telescope Interferometer (VLTI). We derived the size 

5 and the surface thermal properties of (41) Daphne by means of a thermophysical 

6 model (TPM) , which is used for the interpretation of interferometric data for the 

7 first time. Prom our TPM analysis, we derived a volume equivalent diameter for 

8 (41) Daphne of 189 km, using a non-convex 3-D shape model derived from op- 

9 tical lightcurves and adaptive optics images (B. Carry, private communication). 

10 On the other hand, when using the convex shape of Kaasalainen et al. (2002. 

11 Icarus 159, 369-395) in our TPM analysis, the resulting volume equivalent di- 

12 ameter of (41) Daphne is between 194 and 209 km, depending on the surface 

13 roughness. The shape of the asteroid is used as an a priori information in our 

14 TPM analysis. No attempt is made to adjust the shape to the data. Only the 

15 size of the asteroid and its thermal parameters such as, albedo, thermal inertia 

16 and roughness are adjusted to the data. We estimated our model systematic 

17 uncertainty to be of 4% and of 7% on the determination of the asteroid volume 

18 equivalent diameter depending on whether the non-convex or the convex shape 

19 is used, respectively. In terms of thermal properties, we derived a value of the 

20 surface thermal inertia smaller than 50 J s~°'^ and preferably in the 

21 range between and ~ 30 J s~^'^ K~^. Our TPM analysis also shows 

22 that Daphne has a moderate macroscopic surface roughness. 
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1. Introduction 



Information about sizes and shapes of asteroids provides essential cons t raints 



to the history an d formation processes of these bodies (jBottke et al 



2005 



Tanga et al 



20Q9I ). The size distribution of the different subpopulations of as- 



teroids and of the asteroid dynamical families constrain the collisi onal evolution 



proce sses that these bodies have experienced during their history (jBottke et al 



20051 ). Moreover, accurate determination of sizes and shapes is crucial to es- 



timate volumes of asteroids, which allow one to calculate the bulk densities of 
these bodies when their masses are determined by some means. The density and 
internal structure are among the most important characteristics of asteroids, yet 
they are also some of the least constrained. When compared with the densities 
of meteorites - a partial sample of the building blocks of asteroids that survive 
the passage through the Earth's atmosphere - one can deduce the nature of 
asteroid interiors. These physical properties of asteroids reflect the accretional 
and collisional environment of the early solar system. 

The determination of the volumes of asteroids will be particularly important 
in the next future when more asteroid masses are expected to be accurately de- 
rived. For instance, it has been estimated that the masses of slightly more 
than 100 asteroids will be determined to better than 30% (relative accuracy) 
from the gravitational perturbations that these bodies exert on the orbits of 
smaller asteroids thanks to the high accu racy astrometric me asurements of the 



Mouret et al 



20071 ). Yet, the volume 



ESA space mission Gaia (launch in 2013; 
of these bodies are not known with accuracies small enough to allow one to 
calculate meaningful densities. The volumes of asteroids are affected by large 
errors because their true 3-dimensional shapes are generally unknown and thus 
approximated by means of spheres. It can be noted, for instance that by using 
a shape and a spin solution derived f rom lightcurves a nd m utual occultation 



events for the asteroid (22) Kalliope, 



Descamps et al 



(j2008l ) significantly re- 



vised its volume and thus its density compared to previous estimates based on 
a sphere. The error in the estimation of the volume can also be significant when 



large scale topoRraphic co ncavities, known to be present on asteroids (see e.g., 
Thomas et al.l . Il999ll2002l ). are approximated with flat surfaces. 



From the size (D) and the absolute magnitude of an asteroid in the visible 
light (-ff), one can derive the geometric visible albedo (pv) using the formula: 



D(fcm) = 1329 py^^^ IQ-"^^. 



(1) 



The value of the albedo is important to constrain the nature of asteroids: it is 
known, for instance, that asteroids with spectra similar to carbonaceous chon- 
drite meteorites (types CI a nd CM), the so-called C-t ype asteroids, have values 



oipv between 0.03 and 0.10 (jStuart and Binzell . 120041 ): stony (S-type) asteroids, 
rich in silicates such as ol ivine and pyroxene have moderate values of pv (e.g., 
between ^0.15 and ^0.3) ([Stuart and Binzell . 120041 ) . whereas asteroids whose re- 
flectance spectrum is analog to that of ensta tite meteorites are kn own to have in 
general high (> 0.4) albedo values (see e.g. 



Tedesco et al 



1989. and references 



therein) . In general sizes and albedos of asteroids are obtai ned form photometric 
obser vations of these bodies in the thermal infrared (see iHarris and Lagerros , 



2OO2I for a review on the topic) . Models of the surface temperature distribution 



and the cor responding infrared emi s sion are used for the anal ysis of observa- 



tional data ( Delbo and Harris , 



2002 



Harris and Lagerros , 



20021) . In particular, 



thcrmophysical models (TPM) take explicitly into account the effects of thermal 
inertia, spin state, asteroid shape and surface roughness on the calculation of 
asteroids infrared emission. 

One of these parameters, the thermal inertia, a measure of the resistance of 
a material to temperature change, is particularly important. It is defined by 
r = y/pKc, where k is the thermal conductivity, p the density and c the specific 
heat. The valu e of the thermal inertia depends on the material properties (see 



Mueller 



20071 and references therein for a table of thermal inertia values of 
some typical materials) and inform us about the nature of the surface regolith: 
a soil with a very low value of F, for instance in the range between 20 and 50 
J m'^ s~^'^ K~^, is covered with fine dust; an intermediate value (150-700 



J 



-0.5 



K ^) indicates a coar ser, mm- to cm-sized, r egolith as observe d 



on (433) Eros (|Veverka et alJ . l20oiJ u) and (25143) Itokawa (|Yano et al 



20061) 



respectively; solid rock with very little porosit y is known to ha ve thermal inertia 



values of more than 2500 J m ^ s '-'■^ K ^ ([Jakoskv 



19861) . The correlation 



between the value of F and the nature of the soil has been also demonstrated 
from study of the martian surface (see e.g. 



Christensen et al 



20031) . More- 



over, because thermal inertia controls the surface temperature distribution of 
an asteroid, it affects the strength of the Yarkovsky effect. This is the gradual 
drifting of the semi major axis of the orbits of km-sized asteroids caused by the 
asymmetric (with respect to the direction as teroid-sun) emission of the thermal 



94 infrared radiation that carry momentum (see 



Bottke et al 



20061 and references 



therein). This effect plays a role in the delivery o f near-Earth asteroids from 



the main belt ( Morbidelli and Vokrouhlickv 



97 oid families ( Nesvornv and Bottki 



1 



20031 ) ■ in the dispersion of aster- 



20041) ■ and it is a major source of uncer- 



98 tainty in the impact prediction estimation s for potentially hazardous asteroids 



99 (jGiorgini et al 



2002; 



Milani et al 



20091) . Finally, accurate determination of 



thermal inertia is important in the estimation of systematic errors on sizes and 
albedos of asteroids, whe n these parameters a, re determined by means of simple 



thermal models (see e.g 
As shown by 



Spencer et al. 



19891) . 



Delbo et al 



(|2009l ). the Very Large Telescope Interferometer 
(VLTI) of the European Southern Observatory (ESO) can be used to obtain 
measurements of asteroid sizes and shapes. Generally speaking, the VLTI has 
the ability of measuring directly sizes and deriving rough information about the 
shape of asteroids from measurements of the visibility (contrast) of interfero- 
metric fringes. The visibility is a function of the apparent angular extension of 
the body along the projected interferometer baseline. Shape features such as 
large concavities, bilobed shapes and/or presence of satellites, also produce a 
clear signature in the visibility. A sensitive instrument to measure asteroid vis- 
ibilities at the VLTI in the mid-infr ared N-band (8-13 /x m) is the Mid-Infrared 



113 Interferometric Instrument (MIDI; 



Leinert et al 



20031 ). The angular resolv- 
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ing power of the interferometer depends on the length of the basehne. VLTI 
basehnes vary between 16 and 130 m, with theoretical corresponding angular 
resolutions between 130 and 16 mas (milliarcseconds) at 10 ^m. 

Interferometric observations of asteroids with other facilities, such as the 
Fine Guidance Sensors (FGS) of the Hubble Spa ce Telescope (HST, 



Tanga et al 



2001 



Hestroffer et al 



2002 : 



Tanga et al 



20031 ) , demonstrated the capability of 



120 the method of obtaining sizes of asteroids and reconstructing the ellipsoids that 

121 best fit their shape. Given the limiting magnitude of the FGS, only bright 

122 (V<12) and large (^100 km-sized) asteroids were observed by this program. 

123 However, only the large VLTI baselines can overcome these sensitivity and size 

124 limitations by extending the use of the interferometric technique to a large 

125 number of fainter and smaller targets. One particularly interesting feature of the 

126 MIDI instrument is that it also measures the total (non coherent) spectral energy 

127 distribution, /(A), of the source in the 8-13 /im spectral interval. This thermal 

128 infrared data can then be used as a further constraint to derive a steroid sizes , 

129 throug h the application of asteroid thermal models. In their work, iDelbo et al 



( 20091 ) show the first successful interferometric observations of two asteroids 
with MIDI, (234) Barbara and (951) Gaspra. 

In this work, we report on the continuation of the program devoted to mea- 
surement of the physical properties of asteroids from interferometric observa- 
tions in the thermal infrared. In particular, we obtained the first successful 
interferometric observations of asteroids using the Auxiliary Telescopes (ATs) 
of the VLTI. From fluxes and interferometric visibilities measured in the ther- 
mal infrared, we derived the size of the asteroid (41) Daphne, and we studied 
its thermal properties by means of a TPM. 

This work is structured as follows: in section [2] we detail the thermophysical 
model used for the interpretation of MIDI data in terms of asteroid physical 
properties; in section [3] we report the observations and the data reduction pro- 
cess that we adopted; in section [4l we detail the shape models that we used; in 
section [5l we give our results, followed by a discussion in section |6l 
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144 2. Modeling and analysis of MIDI observations 



MIDI is used to coherently combine the infrared light collected by two of 
the four 8-ni Unit Telescopes (UT) or by two of the four 1.8-ni Auxiliary Tele- 
scopes (AT) of the ESO VLTI. The two observables measured by MIDI are the 
photometric flux /(A) and the visibility V{u, v) of the source; where u — B^/X 
and V — By/\ are the spatial frequencies in rad~^ along the x- and y-axis, with 
Bx and By the components, along the two axis, of the interferometer's baseline, 
projected on the plane of the sky. The x and y-axis define the coordinates on 
this plane. We recall that V{u, v) is the Fourier transform of the brightness 
distribution of the source divided by /(A). The visibility can be also indicated 
wit h V^(B/A), where B is a vector of components [B^, By). 



Delboetal 



( 2003) interpreted /(A) and y(B/A) using simple thermal mod- 
els and simple geometric models (disk of uniform brightness and system of two 
disks) in order to derive the size of (951) Gaspra and (234) Barbara. For the 
latter asteroid, a disk of uniform brightness poorly reproduced the observa- 
tions, whereas a b inary model provided a good fit to the data. For this reason 



Delboetal 



(|2009l ) speculated that Barbara has a bilobed shape or a satellite. 
When we have a priori information about the shape and the spin state of an 
asteroid, TPMs can be used to derive the size of the body and to constrain sur- 
face properties, such as albedo and thermal inertia and macroscopic roughness. 
These parameters are explicitly taken into account in the TPM to calculate 
the asteroid's thermal emission, and are adjusted until model fiuxes best fit si- 
multaneously observations obtained at different epochs and wavelengths in the 
thermal infrared. In general, these observations are meas urements of the ob- 
ject's disk integrated therma, l infrared flux /(A) (see e.g., 



2007 



Delbo and Tangal . 120091 ). 



Miiller and Barnes , 



Here, we used a TPM to calculate interferometric visibilities of asteroids in 
the thermal infrared for the first time. Our procedure consists in generating 
images of the thermal infrared emission of the asteroid at different wavelengths 
as viewed by the observer and then obtaining the model visibility and flux for 
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each image: the model flux is calculated by taking the integral of all pixels at 
each wavelength, while the model visibility is calculated as the modulus of the 
Fourier transform of the image, along the projected baseline direction, divided 
by the flux. 

The free parameters of the TPM are adjusted in order to minimize the dis- 
tance between the disk integrated flux /'(A) and visibility V(B/A) of the model, 
and the corresponding observed quantities /(A) and y(B/A). As goodness of 
the fit indicator, we use the reduced x^, namely: 



^ N 



EE 

i=i j=i 



^^fV.iB/X,)-VaB/X,r' 



i=l 3 = 1 



(2) 



where the indexes i and j run over the observation epochs and the discrete 
samples in wavelength at which the visibility and the flux were measured, Nf. 
and Ne are the number of epochs at which flux and visibility are respectively 
measured, Fi and Wi are the number of flux and visibility samples at the i*^ 
epoch, = Y^^=i ^\ I + Ylii=i ^\ V ^^'^ total number of measurements, 
and j and ai- ■ are the uncertainties on the measured visibilities and fluxes, 
respectively. 

The physical parameters of our TPM, are: 

• An a priori information about the shape of the body, described by a mesh 
of planar triangular facets and the spin vector of the asteroid. The shape 
and the spin vecto r are i n general determined by lightcurve inversion (see 



Kaasalainen et al 



20021 for a review), disk-resolved imaging (from e.g.. 



in-situ, the Hubble Sp ace Telescope, and / or ground-based ad aptive optics 



obser vations; see e.g., 



Carrv et al 



2010l) . or radar (see e.g. 



Ostro et al 



20001 ). Our implementation of the TPM allows non-convex shapes to be 



used. Shadowing of facets to the observer due to the body's topography 
is fully taken into account. Mutual heating and light reflection between 
facets within topographic concavities are not modeled in this version of the 
TPM. However, their effect is of second order on the surface temperature 
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determination. The spin vector solution is given by the rotational phase 
(j)Q at a reference epoch to, the ecliptic longitude Ao and latitude /3q of the 
spin axis direction, and the rotation period P. Note that to can be quite 
far in the past. As a consequence, the absolute rotational phase, (f)(t), of 
the asteroid, at a more recent epoch (e.g., the time of VLTI observations), 
can be affected by a significant uncertainty, A^. The latter depends on 
the error crp on P, the value of P, and by how far t is from to- Since 
(f> = (f>o + 27r(t — to)P~"'^, we can write that 



(3) 



neglecting the error on c/jq (which is safe - in general - to assume small). 
When > 10°, then A^ should be treated as a free parameter of the 
TPM. 

• The size of the body. This is described by a factor a that linearly scales 
all vertexes of the mesh. We give the size of the body in terms of the 
diameter of the sphere of equivalent volume D\y = 2 (|^) ^ , where V is the 
volume of the mesh. 



• The bolometric Bond's albedo A. This is related to pv via the relation: 
A = p v{0.29 + 0.68 4 G'), w here G is the slope parameter of the H, G sys- 



tem of 



Bowell et al 



|l989l ). Although the value reported in the MPC for 



(41) Daphne i s 0.10 , the use of the typical default value, namely G=0.15 



(Bowell et al. 



19891) . does not significantly affect our results and we de- 



cided to keep it. 

• The macroscopic surface roughness. This is modeled by adding hemi- 
spherical craters of open ing angle, 7c, and surface density, pc- Following 



Delbo and Tangal (j2009l ) , we used here four preset combinations of 7c and 



Pc spanning the range of possible values of surface roughness. These values 
of macroscopic roughness are given in TablelU inc luding t he cor responding 
value of the mean surface slope, 6, as defined by iHapkd (|l984l ). 
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• The value of thermal inertia, which affects the temperature of each tile of 
the mesh and the temperature distribution inside craters. 

Albedo, thermal inertia, and roughness are assumed constant over the surface 
of the body. 

[Table [T] here] 



In order to compute images of the model thermal emission, the first step is to 
use the TPM to calculate the body's surface temperature distribution, namely 
the temperature of each tile of the mesh and inside craters at any given epoch 
t. This is performed by solvi ng the one-dim ensional heat diffusion equation 



for each tile of the mesh (see 



Mueller 



20071 for more details). The bound- 



ary conditions of the heat diffusion equation are given, at the surface, by the 
input radiative heating from the Sun and irradiation of the heat into space, 
and, at depth, by imposing a zero net heat fiow towards the interior of the 
body. The position of the Sun (the heat source) with respect to each facet is 



244 calculated by applying the inverse coordinates transformations of 



Durech et al 



245 (|2010r) in order to transform the heliocentric position vector of the asteroid into 



the body's mesh reference frame: the transformation is a function of the time t. 
The meth od we adopt to calc ulate the temperature distribution inside craters 



248 is given by 



Emerv et al. 



(119981 ) . We do not explicitly model thermal c onduction 

249 inside craters, and here we use the approximation of iLagerroa (|l998l ). This ap- 

250 proximation cannot be used on the night side. However, this is not a limitation 

251 for the present study because our observations took place at a moderate solar 

252 phase angle (see Table [2]) . Consequently, the fraction of the night side seen by 

253 the observer was negligible. Because of the finite thermal inertia value, the heat 

254 diffusion is not instantaneous and the body temperature distribution depends 

255 on his past illumination history. For this reason, the calculation of the heat dif- 

256 fusion is started ^ 100 rotations (^ a month) before the observation epochs. We 

257 carefully checked that the body temperature distribution stabilizes and is inde- 
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258 pendent of the initial conditions. From the knowledge of the body temperature 

259 distribution, model fluxes are calculated assuming a wavelength independent 

260 emissivity of 0.9, for each tile of the mesh, including craters, in the direction 

261 towards the observer. The emissivity of 0.9 assumed here is typical for the vast 

262 majority of silicate powders. The emissivity of these materi als is within 10 % of 

263 0.9 in the wavelength range between 7 and 14 /im (see e.g.. 



Mueller 



2007, and 



264 references therein). Then the three-dimensional mesh is projected on the plane 

265 of the sky to create a two dimensional image of the asteroid. In order to have 

266 a reasonable image size, while remaining close to the maximum resolution of 

267 the MIDI instrument, we sampled the images with a resolution of 4 mas/pixel 

268 ( 5 mas/pixels for the convex shape of 41 Daphne; see section [S]). Our results 

269 are robust with respect to changes of the pixel scale. The value of each pixel is 

270 calculated from the flux of the facet of the mesh on which the pixel is projected 

271 to: namely it is the facet flux multiplied by the area of the pixel and divided 

272 by the projected area of the facet. Fig. [T] shows two images of (41) Daphne 

273 obtained along the lines described above, from the two shape models described 

274 later. These images are created for each wavelength. 

275 

276 [Fig. [T]here] 

277 

278 Then the integrated flux and visibility of the model are given by: 

/'(A,) = a2^0'(x,y,A,) (4) 



y'lB\^FT{0'{^,y-,X,))_a&i^) 



A,7 ^'(A,) /'(A,) 



-i2T7 



= 77^)T.[T.O'i-'^y'^^.)je^' (5) 

279 where FT is the Fourier Transform operator applied to the brightness distribu- 

280 tion of the asteroid O', (x', y') is a system of coordinates rotated by an angle PA 

281 with respect to the {x, y) frame. The y' axis thus coincides with the direction 
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282 of the baseline (see Fig. [5]) so that we replace the baseline vector, noted B, by 

283 its modulus, noted _B, in the expression of the model visibihty above. Nyi is 

284 the number of pixels in the image along the y' axis. In Eqs. |3]and[Sl we have 

285 explicitly written the dependence of the flux and the visibility on the mesh scale 

286 factor a. It is important to note that changing the scale of an object by a linear 

287 factor a is equivalent to multiply the baseline length by a factor a, or inversely. 

288 

289 [Fig. H here] 

290 

291 The free parameters of the TPM are: size (which is varied through the scale 

292 factor a), thermal inertia, and macroscopic surface roughness. The fit procedure 

293 involves calculation of /j'(Aj) and Vl{B/\j) at each observing epoch and each 

294 wavelength Aj, for a number of discrete values of F and macroscopic surface 

295 roughness. Note that the dependence of the on a is trivial, as it can be seen 

296 from Eqs. |3]and[S] Thus, the best value of a for each discrete value of F and 

297 each roughness model can be found by minimizing the Eq. [2] Then, the location 

298 of the minimum as function of F gives the best-fit asteroid surface thermal 

299 inertia for each roughness model. Eventually, the value of a at F-niinimum is 

300 used to determine the best- fit value of Dy . 

301 In some cases, the correction to the rotational phase has to be treated as 

302 a free parameter of the TPM. In this case, the location of the minimum ^f^, as a 

303 function of for different roughness and thermal inertia values, gives the set 

304 of best-fit physical parameters for the asteroid. An optical lightcurve obtained 

305 quasi-simultaneously could help reduce A^. 

306 

307 In the next section, we describe our observations of the asteroid (41) Daphne 

308 and the related data reduction. 
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309 3. Observations and data reduction 



310 The observations of (41) Daphne were carried out in service mode, on 2008 

311 March 12 and 14 . We adopted the typical observing s equence with MIDI, 



which is extensively described by 



Przygodda et al 



(120031) . Using two ATs in 



the EO-GO baseline configuration {B = 16 m), we acquired four visibility obser- 
vations for (41) Daphne, two on March 12 and two on March 14. Table [2] reports 
the observational circumstances. The telescopes and the delay lines of the inter- 
ferometer were tracked at the rates predicted from the ephemerides of the target. 

[Table. [2] here] 



320 Our observations included a mid-infrared photometric and interferometric 

321 calibrator chosen from the ESO database, namely IID123139. The abs olutely 



Cohen 



322 calibr ated infrared spectrum of the calibration star was taken from 

323 (|l999f ). Extraction and calibration of the visibility mea s ureme nts of (41) Daphne 

Delbo et all |2009l ). The fiux and vis- 



324 were performed using the same method as[ 

325 ibility measurements of (41) Daphne are shown in Fig[3](or identically in Fig|4]). 

326 The estimation of error bars constitutes a difficult issue when reducing MIDI 

327 data. In the most common case, when the 'high-sensitivity' mode is used, the 

328 photometry is acquired about 3 to 6 minutes after the fringes are recorded. 

329 Therefore, the measured value of /(A) does not correspond to the flux of the 

330 source at the time of th e fringes r e cordi ng. This leads to typical uncertainties 

331 of about 10 to 15% (see 



Chesneaul . 



20071 for more details). This error strongly 



332 depends on the atmospheric conditions during the night. An estimation of the 

333 error bars can thus be obtained by computing the RMS of the visibilities val- 

334 ues of several calibrators observed closely to the source (or identically several 

335 observations of the same calibrator). We used this procedure for the two first 

336 visibilities, and considered a typical relative uncertainty of 15% for the two 

337 others since only one calibrator was observed during the second night. 

338 [Fig. H and Fig. H here] 
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340 In addition to the interferometric observations, we obtained an optical lightcurve 

341 of (41) Daphne, in order to better constrain the rotational phase at the epoch 

342 of the VLTI observations. In October 2009, during the tests of the 0.4 m 

343 telescope for the An tarctic Search for Transiting ExoPlanets (ASTEP) project 



344 ( Daban et al 



20101) at the Observatoire de la Cote d'Azur in Nice, France, we 



345 performed time-resolved CCD photometric observations of this asteroid in the 

346 visible. Our observations spanned a period of about 3 h. We used the technique 

347 of differential photometry, which makes use of stars present in the same CCD 

348 frame as the target, to allow for accurate removal of systematic effects such as 

349 atmospheric extinction variability. The telescope was tracked at the sidereal 

350 rate. The reduction of the CCD frames consisted in the conventional bias re- 

351 moval and flat-fielding, which was performed by using bias and flat calibration 

352 frames. The corresponding lightcurve is shown in Fig. [51 where the times of the 

353 observations were light-travel subtracted. 

354 

355 [Fig. [5] here] 

356 

357 In the next section we describe the different shape models used for (41) 

358 Daphne, and how the TPM analysis was performed using this information). 



4. Shape models and TPM analysis 



360 For sake of comparison and evaluation of the results variability, two shape 

361 models were used in this study: 

362 1. a convex mesh was downloaded (March 2010) from th e Database of Aster- 



oid Models from Inversion Techniques or DAMIT (see 



Durech et al 



2010 



andhttp: //astro . troja.mf f . cuni . cz/projects/asteroids3D/web .php?page=project_main_page) 

This is a database of three-dimensional models of asteroids derived by 
solving the inverse problem of determining the object's shape and its ro- 
tational state from optical lightcurves. It also lists some models derived 
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from a combination of adaptive optics images and lightcurves. The convex 
shape of (41) Daphne contains 1022 vertexes and 2040 triangular planar 
facets. The pole solution is also given: Ao=198'' /3o=-32'', P=5.98798 h, 
00=0° at the reference epoch to=2444771. 79382 (JD); 
2. the second shape model of (41) Daphne, provided by B. Carry (private 
communicati on]^, i s non-convex and was derived from the KOALA method 
2010; 



( Carrv et al 



Kaasalainen . 



2011 ). We recall here that this method 



375 allows one to derive the size and the shape of an asteroid from the inversion 

376 of data obtained by different observational techniques: namely, photomet- 

377 ric lightcurves and adaptive optics images in the case of (41) Daphne. The 

378 corresponding pole solution is: Ao=198° /?o=-31°, P=5. 987980 h, ?!)o=0° 

379 at the reference epoch to= 2444771.79382 (JD). 

380 Note that, while the size of the first shape model is arbitrary, the non-convex 

381 shape, derived from the KOALA method, has a volume equivalent diameter of 

382 D\j = 185 ± 5 km. We call it 'the nominal size' in the following. The pole 

383 solution of the convex model of (41) Daphne was obtained from inversion of a 

384 set of 23 optical lightcurves obtained in the period 1976 - 1988. Assuming = 

385 December 21, 1988 and crp=0. 00001 hours, we thus find that A0 17° at the 

386 epoch of our VLTI observations, given the rotation period of 5.98798 hours for 

387 this object. As a consequence, we attempted to better constrain the rotational 

388 phase of (41) Daphne using our recently obtained lightcurve shown in Fig. [5l 

389 To do so, we calculated synthetic visible lightcurves using the convex shape of 

390 ( 41) Daphne and attempted to fit the model to the observed lightcurve using 

391 the value of as a free parameter. A geometric scattering model was adopted: 

392 namely the visible reflected flux from each facet is proportional to the projected 

393 area of the facet, and it is required that the facet is also illuminated by the 

394 sun. Fig. [5] also shows the best fit model lightcurve that we obtained using the 

395 convex model for a correction value of — (—0.2 ± 0.3)°. The correction to 



^ See also his PHD thesis manuscript available at [http : //benolt . carry . free . f rZ] for fur- 
ther information. 
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396 the nominal rotational phase is thus negligible. 

397 For the non-convex model, given that the KOALA model uses recent observa- 

398 tions from 2008 March 28, and that the period is more accurate (by an order 

399 of magnitude), the error on the rotational phase given by the pole solution is 

400 negligible as well. As a consequence, for the two shape models, the rotational 

401 phase was not treated as a free parameter. The TPM was run for each roughness 

402 model (see Tabled]), and each value of thermal inertia, namely 0, 5, 10, 25, 50, 

403 75, 100, 150, 200, 300, 400 and 500 J m^^ 5-0.5 ^hen the fit procedure 

404 described in section [2] was applied to the measured fluxes and visibilities. 

405 In the next section we describe and discuss the results obtained from the 

406 application of the TPM to the observed visibilities and fluxes of (41) Daphne. 

407 5. Results 

408 5.1. Convex shape model 

409 Fig. |6] shows our best-fit estimator as a function of F for the four different 

410 roughness models, in the case of the convex shape. 

411 

412 [Fig. [6] here] 

413 

414 We note that a surface with a low or no macroscopic roughness and a value of 

415 thermal inertia < 100 J m^^ s^'^'^ K^^ gives the best fit to the observations. 

416 In particular the minima of the 'no roughness' and 'low roughness' models are 

417 at F=48 and 8 J m^^ s^"-^ K~^, respectively. The corresponding values of I?v 

418 are 209 and 194 km, with associated pv values of 0.057 and 0.067, respectively. 

419 In Fig. 13] we plot the visibility and flux of the best- fit model ('no roughness', F = 

420 48 J m^^ s^^-'^ K^^, Dy — 209 km, pv — 0.057), in addition to the measured 

421 fluxes and visibilities of (41) Daphne. We note that our model represents well 

422 the observed flux except for the third observing epoch, where the fiux of the 

423 TPM is greater than the measured one (by roughly 20% at 13 /im). This 'offset- 

424 like' mismatch may likely come from an underestimation of the source flux by 
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425 MID I, especial l y arou nd 13 /im where the signal to noise ratio is generally quite 



426 low ( Chesneau 



20071 ). The photometry measurement of the source may thus be 



427 degraded by a bad estimation and suppression of the thermal bac kground (and 

428 its flu ctuations), which is dominant in the mid- infrared (see e.g., 



Perrin et al 



429 l2003f) . Moreover, the underestimation of the flux at the third epoch implies a 

430 similar offset-like mismatch between the corresponding measured and calculated 

431 visibilities at the same epoch. 

432 In order to estimate the statistical uncertainty affecting the fit parameters 

433 r and D\y, a Monte-Carlo analysis was performed. To this end, 200 normally 

434 distributed flux and visibility values per observation were generated at each 

435 wavelength, with average and standard deviation matching the measured fluxes 

436 and visibilities within their respective 1 — uncertainty. For each set of fluxes 

437 and visibilities, a fit of the model, as described in section [21 was performed. 

438 Then, we took the standard deviation of the F and D\/ values at the minimum 

439 as the 1 — cr uncertainty on our best fit values of thermal inertia and volume 

440 equivalent diameter. As a result we find F = 48±5 J to~^ s~" ''^ K"^ and _Dv = 

441 209±1 km as the best fit solution for a model using the convex shape without 

442 roughness; and F = 8±8 J s^"-^ K^^ and Dsj = 194±2 km as the best fit 

443 solution for a model using the convex shape with a low roughness. 

444 5.2. Non-convex shape model 

445 Figure [7] shows our best-fit estimator x'^ of the TPM, calculated in the case 

446 of the non-convex shape. 

447 

448 [Fig. [7] here] 

449 

450 In this model with a low or medium roughness and a thermal inertia 

451 < 100 J m^^ s^'^'^ gives the best fit to the observations. In particular, the 

452 minima of the 'low roughness' and 'medium roughness' models are at F— 9 and 

453 J m^^ s^°-^ K~^, respectively. The corresponding value of are 189 km and 

454 182 km, respectively, with pv~ 0.07. The solution with F=0 J s^'^-^ K^^ is 
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455 non physical and thus we reject it. Although 100 km-sized Main Belt aster- 

456 oids are known to have low thermal i nertia, their values are in ge neral larger 



Mueller and Lagerros , 



19981 and refer- 



457 than 10 J m^^ s^°-^ (see e.g., 

458 ences therein). We also used the Monte-Carlo procedure described above to 

459 estimate the uncertainties on the fit parameters. These uncertainties result to 

460 be of 1 J m^^ s^'^-^ K^^ for the thermal inertia and 1 km for the volume 

461 equivalent diameter. We show in Fig. 0] the measured fluxes and visibilities of 

462 (41) Daphne along with the best fit of the non-convex case ('low roughness', 

463 r = 9 J j^-i^ ^ ;^39 i^j^^ ^ 0.070). Our model, when used 

464 with the non-convex shape, also provides a good match to the observed flux. 

465 We note that the flux prediction at the third epoch is as well higher than the 

466 measured one. This effect may also be due to an underestimation of the source 

467 flux by MIDI, as described above. Moreover, the non-convex model seems to 

468 better flt the measured visibilities than the convex one. This result is consistent 

469 with the fact that the non-convex shape model corresponds to a more realis- 

470 tic and detailed representation of (41) Daphne, as previously indicated by B. 

471 Carry (private communication)!^ who better reproduced an occultation proflle 

472 of (41) Daphne with this shape model. We can note however a slight mismatch 

473 between the model predictions and the measurements at the fourth epoch. The 

474 mismatch is located around ~ 6 as~^ (or A ~ 9.6/im) i.e. close to the atmo- 

475 spheric ozone absorption feature. This time-dependent absorption feature can 

476 cause fluctuations and lower the signal to noise ratio on both the coherent flux 

477 and the photometry. Due to its relative variability, this feature is frequently im- 

478 perfectly removed during th e reduction process of interferometric observations 



479 (see e.g.. lTubbs et al.l . I2004I ). All the results are summarized in Table [3l 



[Table H here] 
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483 6. Discussion 



484 6.1. Size and geometric visible albedo 

485 The best-fit value of Dsj obtained from our TPM analysis of MIDI data us- 

486 ing the non-convex shape model, namely 189 ± 1 km, presents a discrepancy 

487 with the nominal value (185 ± 5 km) of about 4 km. This is about 4 times 

488 greater than the statistical uncertainties implied by the photometric and vis- 

489 ibilities measurements and estimated via our Monte Carlo method. However, 

490 this discrepancy represents a relative accuracy of 2% on the diameter that ap- 

491 pears remarkably good for a typical size determination of asteroids. If the two 

492 equally plausible solutions in terms of our best-fit indicator, i.e. low roughness 

493 - r=9 J m^^ s~°-^ K^^ and medium roughness - r=0 J m^^ s^"-^ K~^, are 

494 considered in order to estimate the model error on the value of D\j, we found 

495 that the latter is 7 km i.e. 4% in relative accuracy. 

496 When the convex model is used as an a priori shape for the TPM, the value 

497 of Dsj resulting from the analysis of MIDI data is between 194 and 209 km (with 

498 a mean value of about 200 km), depending on whether a low-roughness or a no- 

499 roughness model, respectively, is assumed for the surface of Daphne. Because 

500 both solutions are equally good in terms of our best-fit estimator, it is clear that 

501 the systematic model uncertainty dominates the final size determination. This 

502 systematic relative uncertainty is of the order of 7% on the diameter. 

503 It is also evident that the best fit values of the equivalent volume diameter 

504 obtained when using the convex shape are larger than the nominal value of 

505 Dy of the non convex shape. This is because the latter has a smaller volume 

506 (due to the presence of topographic concavities) for a similar projected size on 

507 the plane of the sky. This volume overestimation when the convex shape is 

508 used, clearly indicates that the imp osed condition of convexi t y, imp osed by the 



20011 ). introduces 



509 lightcurve inversion technique (see iKaasalainen and Torppa , 

510 a significant systematic bias on the size when large concavities are present on 

511 the asteroid surface. 

512 The geometric visible albedo values derived from the TPM with the convex 
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513 and non-convex shapes are in the range between 0.05 and 0.07. This is in 

514 agreement with the C-type taxonomic classification of (41) Daphne. 



515 6.2. Thermal properties 

516 We obtained for the first time good constraints on the determination of 

517 macroscopic roughness for both shape models, although the results slightly dif- 

518 fer. While a very low roughness is preferred in the convex shape low 

519 or moderate roughness gives the best-fit in the case of the non-convex one. We 

520 estimate that the corresponding mean surface slope values should lie between 

521 10 and 29° for (41) Daphne. A high macroscopic roughness is discarded, which 

522 is quite surprising given the important surface roughness expected for large as- 

523 teroids. Large main-belt asteroids such as (1) Ce r es and (2) Pallas are expected 

524 to have rough surfaces as found bv lSpencer et al.l (|l989f ) from the analysis of the 

525 thermal infrared emission of these bodies. It is also known that the observed 

526 zero-phase thermal emission of the Moon is well r eproduced by a rough-surface 

527 model, i.e. pc = 0.64 and 7c = 90° (|Spenceii . ll990l ). The possibility to constrain 

528 the macroscopic roughness in the thermal infrared is very interesting because 

529 it allows, in principle, a more accurate determination of the asteroids thermal 

530 inertia from TPM modeling. 

531 Our TPM analysis indicates that (41) Daphne has a thermal inertia value 

532 certainly smaller than 50 J m^^ s^^'^ K^^, and likely in the range between 

533 1 and ^ 30 J s^°-^ K^^, indicating that a layer of very fine regolith 

534 covers its surface. This is in agreement with the thermal inertia values already 

535 measured on main-belt asteroids lar ger than 100 km in diameter . Indeed, from 



536 mid- and far-infrared observations, 



Mueller and LagerrosI (jl998l) derived very 



low thermal inertia values, between 5 and 25 J m ^ s K ^, for some large 
Main Belt asteroids includin g (1) Ceres, (2) Pallas, (3) Ju no, (4) Vesta, and 



539 (532) Herculina. Moreover, 

540 inertia of about 15 J s 



Miiller and Blommaerti (|2004l ) derived a thermal 



-0.5 



K ^for (65) Cybele. More recently, from the 



541 direct measur e ment of the cooling during shadowing events in a binary system, 
Mueller et al.l (|2010l ) measured a thermal inertia of 20± 15 J s~°-^ K^^ioi 
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the large binary Trojan {D « 100 km), (617) Patroclus. 

Several works were recently devoted to the determination of the thermal 
inertia value of the asteroic 
mal i nfrared data (see e.g., 



Delbo et al. 


2007 


Delbo and Tanea 


2009 


Mueller 



20071 and references therein). In general, these works clearly showed that a 
good fit of the model to the observations can be obtained by using any rough- 
ness model. Because the determination of the value of the thermal inertia is a 
function of the roughness model adopted, any constraint on the latter parameter 
will allow improving the determination of thermal inertia. In order to show the 
effect of the visibility on the determination of a simultaneous solution for both 
thermal inertia and surface roughness, we also performed a TPM analysis of 
the disk integrated thermal infrared data of (41) Daphne, using the non-convex 
shape and neglecting the visibility measurements. 



[Fig. H here] 



As it can be seen from Fig. [51 any roughness model and any value of the 
thermal inertia give a value of smaller than 1, which implies that any value 
of these parameters allows the TPM to fit equally well the infrared flux. In 
order to reduce this degeneracy, it is known that thermal infrared photometric 
observations are needed at different illuminations and viewing geometries in 
order to constrain the asteroids thermal properties from a TPM modeling of 
photometric observations only. 



567 7. Conclusion 

568 We have obtained the first successful intcrferometric observations of asteroids 

569 in the thermal infrared using the ATs of the ESO VLTI. We observed the asteroid 

570 ( 41) Daphne using the MIDI instrument and the 16ni-long baseline EO-GO. 

571 We developed a thermophysical model (TPM) for the analysis of interfero- 
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572 metric observations of asteroids in the thermal infrared, with the aim of deriving 

573 information about size and thermal properties. 

574 We applied our TPM to the MIDI observations of (41) Daphne: our results 

575 indicate that Daphne has a volume equivalent diameter between 194 and 209 

576 km, depending on the surface roughness, if a convex shape model derived from 

577 lightcurve inversion is used as an a priori constraint on the shape of the asteroid. 

578 Since the nominal size, attached to the non-convex KOALA model, is 185 ±5 km 
57. (B. Carry, private communicationjSl, our results confirm that the assumption of 

580 convexity introduces a systematic bias on the size determination when important 

581 concavities are present on the asteroid's surface. In contrast, if the non-convex 

582 shape model is used, the volume equivalent diameter obtained from the TPM is 

583 189±1 km, i.e. very close to the nominal value of 185 ± 5 km. We estimated our 

584 model systematic uncertainty to be of 4% and of 7% on the determination of 

585 the asteroid volume equivalent diameter depending on whether the non-convex 

586 or the convex shape is used, respectively. 

587 Our TPM analysis also showed that the macroscopic surface roughness can 

588 be constrained by interfcrometry, thanks to the angular resolving power offered 

589 by the VLTI and which allows to resolve the temperature distribution on the 

590 asteroid surface. In particular, using both shape models of (41) Daphne, we 

591 found a moderate to low roughness (see Table [T|); 'high roughness' models are 

592 discarded by our analysis. With such a constraint on the macroscopic roughness, 

593 the TPM results indicate a very low thermal inertia for (41) Daphne, certainly 

594 smaller than 50 J m^^ s^"'^ K^^. This confirmed previous results indicating 

595 that the surface of asteroids with sizes larger than 100 km has a low thermal 

596 inertia. As shown by this work, the possibility of constraining the macroscopic 

597 roughness is important in the prospect of an accurate modelling of the thermal 

598 infr ared emission of asteroids, and especiall y thermal inertia determination (see 



also 



Muelleil . 



2007 



Delbo and Tanga . 



20091) ■ 
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Tables and Table Captions 



Roughness model 7c Pc 



no roughness 


0° 


0.0 


0° 


low roughness 


45° 


0.5 


10° 


medium roughness 


68° 


0.8 


29° 


high roughness 


90° 


1.0 


58° 



Table 1: The four roughness models used in the application of the TPM to the MIDI data; 7, 
and pc respectively correspond to the crater opening angle and the crater density, w hile 9 is 
the co rresponding mean surface slope accord ing to the parameterization introduced by 



1 1984 ) (see text and also iDelbo et al.L 120071 , for further details) 



Hapke 



Date 


UT 


r (AU) 


A (AU) 


a (deg) 


B (m) 


PA (deg) 


Tag 


2008-03-12 


05:19:14 


2.0866 


1.1920 


15.8 


13.6 


66.6 


Di 


2008-03-12 


06:24:35 


2.0864 


1.1917 


15.8 


15.4 


70.6 


D2 


2008-03-14 


04:19:27 


2.0834 


1.1769 


15.0 


11.6 


61.6 


D3 


2008-03-14 


04:32:48 


2.0834 


1.1769 


15.0 


12.2 


63.4 


D4 



Table 2: Observational circumstances and interferometric parameters of the (41) Daphne 
observations, r and A are the heliocentric and geocentric distances, respectively, while a is 
the solar phase angle. B and PA are respectively the length and the position angle of the 
baseline projected on sky. The last column contains a tag associated with each observation. 
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Shape Roughness model ^ V Dv vv 

(J s-0.5 ^-1) (kjn) 

Convex No roughness 3.5 ± 0.3 48 ± 6 209 ± 1 0.057 ±0.009 

Low roughness 3.7 ±0.3 Stg" 194 ± 2 0.067 ±0.011 

Non-convex Low roughness 3.4 ±0.3 9±1 189 ± 1 0.070 ±0.011 

Table 3: Results of the determination of physical properties of the asteroid (41) Daphne, using 
the TPM. The is our best-fit estimator as described by Eq[2l F is the thermal inertia, Dv is 
the spherical volume equivalent diameter, and py is the geometric visible albedo. The errors 
are within l-cr. 
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Figure captions and figures 
Fig [1] caption : 

Image of the asteroid (41) Daphne created from the TPM, using the convex 
shape model (left image) and the non-convex one (right image) . The gray level 
is proportional to the emitted thermal infrared flux. 
Fig [2] caption : 

Illustration of the geometric parameters involved in the calculation of the syn- 
thetic visibility and flux from a TPM image. 
Fig [3] caption : 

Left panels: measured thermal infrared fluxes (with error bars) between 8 and 
13 /im of (41) Daphne, and the corresponding best-fit synthetic infrared fluxes 
(solid lines) derived from the TPM in the case of the convex shape; right panels: 
measured interferometric visibilities plotted, in N band, as a function of angular 
frequency, and the corresponding synthetic visibilities of the TPM (solid lines). 
The best-fit model represented here is : 'no roughness', F = 48 J to~^ s~°-^ K~^. 
The tags Di, D2, D3 and D4 indicate the observing epoch in the chronological 
order (see Table [2]). 
Fig [4] caption : 

As for Fig. 131 but the best-fit model represented here, in the case of the non- 
convex shape, is : 'low roughness', F = 9 J m^^ s^"-^ K^^. The tags Di, D2, 
Z?3 and D4 indicate the observing epoch in the chronological order (see Table [2]). 
Fig [5] caption : 

Diamonds: lightcurve obtained from CCD photometric observations in the vis- 
ible of (41) Daphne in October 2009 during the tests of the ASTEP 0.4 m 
telescope at the Observatoire de la Cote d'Azur, Nice, France. These observa- 
tions spanned a period of about 3 hours. Solid curve: corresponding best-fit 
lightcurve from the convex model. 
Fig [6] caption : 

Plot of (see Eq.[2]), calculated from the TPM in the case of the convex shape, 
as a function of thermal inertia F, for the four roughness models (see Table [1]). 



35 



Fig [7] caption : 

Plot of (see Eq. [2|) , calculated from the TPM in the case of the non-convex 
shape, as a function of thermal inertia F, for the four roughness models (see 
Table [I]). 

Fig \8\ caption : 

Plot of (see Eq. [5]) calculated from the TPM, using only the flux measure- 
ments of (41) Daphne. This is represented in the case of the non-convex shape, 
as a function of thermal inertia F, for the four roughness models (see Tabled]). 
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Figure 1: 
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Figure 4: 
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